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Within the framework of Lorentz-violating extended electrodynamics, the Dirac equation for a 
bound electron in an external electromagnetic field is considered assuming the interaction with 
a CPT-odd axial vector background 6 M . The quasi-relativistic Hamiltonian is obtained using a 
1/c-series expansion. Relativistic Dirac eigenstates in a spherically-symmetric potential are found 
accurate up to the second order in bo. bo-induced CPT-odd corrections to the electromagnetic dipole 
moment operators of a bound electron are calculated that contribute to the anapole moment of the 
atomic orbital and may cause a specific asymmetry of the angular distribution of the radiation of a 
hydrogen atom. 
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I. INTRODUCTION 

O; 

The Standard Model is currently proved with a convincingly wide set of experiments. Nevertheless, it essentially does 
not include a quantum description of gravitation. The quantization methods adopted in the Standard Model do not 
allow a self-consistent quantization of General Relativity since the subsequent theory occurs to be nonrenormalizable. 
Thus the essence of the theory accounting for the effects taking place at the Planck scale of energies (Ep\ = 10 19 GeV) 
where quantum gravity plays a major role, still remains obscure. At the same time, there exist some candidates for 
such Fundamental theory, string theory for instance, taking the form of the Standard Model in the low-energy limit. 
Planck energies being far from experimental attainment, the Standard Model Extension (SME) was elaborated. It 
i-C . is an effective theory (applicable at the energies E <C Ep\) formulated axiomatically as a set of corrections to the 
Lagrangian of the Standard Model fulfilling some 'natural' requirements [l], Q such as observer Lorentz invariance, 
, 4-momentum conservation, unitarity, and microcausality. In what follows, we will focus on a subset of the SME 
' referred to as the minimal SME in flat Minkowsky spacetime that also requires local SU(3)c x SU(2)j x f7(l)y 
gauge invariance and power-counting renormalizability. A spectacular feature of such requirements is that they 
^ ■ reduce the diversity of possible corrections down to a finite number of them. Each correction term consists of a 
complex (pseudo)tensor constant (SME coefficient) contracted with conventional Standard Model fields and their 
spacetime derivatives. These constants are believed to stand for vacuum expectation values of the fields featuring in 



the hypothetic Lorentz-covariant Fundamental theory and condensed at low energies due to the spontaneous symmetry 



o 
o 



X 

5-H 



breaking mechanism. Indeed, it has been shown recently that such Lorentz symmetry breaking can occur in some 
theories beyond the Standard Model [1, 0, d, @ leading subsequently to the SME. The SME can thus be used to 
reduce the complexity of these theories and related calculations in the low-energy limit. It also provides a standard 
for representation of data obtained in experiments searching for Lorentz violation. 

Recently, a number of theoretical researches have been performed aiming at investigating the vacuum structure 
of this model (see, e.g. 0, B S Ell EI), an d to study the assumed violation on various high-energy processes 
[HI, [IH, [li], [l5[ . This search also seems quite promising in atomic physics [ljl [l7| ■ For instance, specific types of 
Lorentz violation may cause spatial parity violation in electrodynamics at tree level. P-parity violation effects in atomic 
systems within the conventional Standard Model have been thoroughly studied either theoretically or experimentally 
in the past four decades [HI, [h| [2(| HH, [HJ . Such effects are caused by weak interaction and include resonant dichroism 
of atomic gas, permission of conventionally forbidden quantum transitions etc. Much the same effects are expected 
within the SME. 

Until today however, studies of atom within the SME have included only spectroscopic predictions using the 
perturbation theory with respect to SME coefficients [l6|, [l?], HH 24] . Direct solution of the atomic eigenstate problem 
would make it possible to study radiative properties of the atom, too. This paper is devoted to an analysis of dynamics 
of a charged fermion in an external electromagnetic field within extended electrodynamics with a background axial 
vector b^. Approximate methods are used to implicitly solve the eigenstate problem in a central electric field; for the 
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Coulomb field, an explicit solution is obtained (see section FV]). The quasi-relativistic approach is also employed to 
obtain the corrections the conventional Schrodinger-Pauli-Dirac equation acquires in the background of bg (sections 
PH IIV[) . Moreover, in view of recent publications [IS US H3| j some other types of couplings and the corresponding 
fcrmion eigenstates are also discussed. As an example of applying these results, an effect of fro-induced asymmetry of 
the angular distribution of spontaneous radiation of a polarized hydrogen atom is demonstrated (section IVI[) . 

In addition, a polarized hydrogen atom is shown to have a nonzero anapole moment originally introduced in [l8j |. 
This characteristic is specific for parity-nonconserving systems interacting with electromagnetic field. For example, 
due to weak interaction loops, neutrinos can possess such a moment, and it is the only electromagnetic characteristic 
that maintains for Majorana neutrinos [28[. 



II. THE MODEL 



We will restrict our consideration to a specific case of extended electrodynamics of electrons and photons within 
the SME (further refereed to as extended QED) with the Lagrangian 

C = -\f^ + i, [^f*D„ -m e - b^l^j i>, (2.1) 

D M = d^ + ieA^x), 75 = -i 7 ° 7 V7 a , (2.2) 

where e, m e are the electron charge and mass, respectively; = j et, r}, and 6 M is a constant CPT-odd axial vector. 
Present constraints on for electron are the following [ljl, H7L l29f : 

I bo I < l(T 2 eV, (2.3) 
\b\ < 1(T 19 eV, (2.4) 

while the constraints for nucleons are some orders more stringent. We will use one-particle approximation in the 
framework of relativistic quantum mechanics. Recalling (|2.ip . one can obtain the Hamiltonian for an electron in an 
external electromagnetic field 

H{t) =a P + j3m e + eA (r, t) + 6075 + b ■ S, (2.5) 

where P = p — eA(r, i) and a. = 7°7, j3 = 7 , and S = —0:75. In Dirac spinor representation, spatial parity operator 
takes the form: 

P£(r,<) = 7°£(-r,i), (2.6) 

hence P f = P, P^P = P 2 = 1. Hamiltonian ([2~5|l commutes with P if A (r,t) = A (-r,t), A(r,t) = -A(-r,t), 
and bo = 0. In particular, the presence of b can violate the P-parity of Hamiltonian (|2.5[) in a spherically-symmetric 
field A^ = {4>(r), 0}, in the Coulomb field of an infinitely heavy nucleus for example, with 

<t>{r) = ~. (2.7) 
47rr 

Charge conjugation of Hamiltonian (|2.5p only changes the sign of electric charge e to the opposite; so hydrogen and 
anti- hydrogen atoms possess equivalent dynamics even if b^ ^ 0. Due to these facts our primary interest concerns the 
investigation of unusual properties of a hydrogen atom induced by the presence of a nonzero constant 6q. 



III. l/c 2 -APPROXIMATION FOR THE DIRAC EQUATION IN THE fe M ^ CASE 

The quasi-relativistic approximation assumes an expansion into a series with respect to 1/c. It is thus necessary in 
this section to rewrite the Hamiltonian (|2.5[) in the CCS system of units with the speed of light c 7^ 1 

H = ca ■ P + 0m e c 2 + eA + cb tl5 + b ■ S, (3.1) 

where P = p — |A(f,i) and bt = bo/c so that bt has dimensionality of momentum. Consider the extended Dirac 
equation with the Hamiltonian (|3.ip in a non-stationary external field A^{x): 

ifi ^M = H{tMr,t), (3.2) 

J ^{r,t)iP(r,t)d 3 r = 1. (3.3) 
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Following the standard method (see, e.g. [20, [31(), let us shift the energy by means of a unitary transformation 



,m e c 2 , I lu 



^^"V J' (3 - 4) 



In terms of 2-component spinors u and v, the Dirac equation (|3.2|) takes the form 



A cA \ fu 
cA A - 2m e c 2 ) \v 



0, (3.5) 



d 

A = <t • P - b t , \ = eA + a-b-ih—, (3.6) 

at 

where <x denotes a vector of the three Pauli matrices, and the Dirac matrices are taken in the standard representation: 

Consider an electron in a state with a positive sign of energy (this in fact does not imply that the electron possesses 
a definite energy). External fields are assumed to be weak enough and to have frequencies much smaller than m e c 2 /H 
so that E, H, P, ih4r — eAo = O(c ), when acting upon u, v. 



In this section, we also assume that 



b t = b /c=O(c°), b = 0(c°). (3.8) 



Contrary to the conventional electrodynamics, within the context of which the Gaussian units are usually used, the 
order of 6 M in 1/c is quite ambiguous. In electrodynamics, certain powers of 1/c can be assigned to the fields E, H, 
which result in a hierarchy of electromagnetic effects having different orders in 1/c. For instance, radiative processes 
are at least of the third order in 1/c, hence the l/c 2 -approximation is worth considering. In contrast, the physical 
origin of b^ is not yet finally established, and hence, we use the convention (|3.8p . due to the symmetry between er • P 
and bt, both entering (|3.5p . Together with (|3.6[) , this convention implies 

A,A = 0(c ), (3.9) 

when acting upon u, v. Then the second line of (|3 . 5[) gives: 



' l + TT^Tl Au + 0(1/c 4 ), (3.10) 



2m e c \ 2m e c 2 

and v is thus suppressed, compared with u, for nonrelativistic positive-energy solutions. On the other hand, the 
square of the norm 



4>^d 3 r = 1 = 0(c°), (3.11) 
= / ( u tu + v*v)dPr = { u\l + 0{l/c 2 ))ud 3 r, (3.12) 



consequently, u = 0(c ) and, due to (|3.10[) . v — 0(l/c). This, in addition, results in suppression of the terms in 
the quasi-relativistic Hamiltonian stemming from the block-off-diagonal part of the matrix in (|3.5[) . in particular, the 
terms containing bo. As a result, the quasi-relativistic Hamiltonian will contain leading-order 6o-induced contributions 
proportional to b t = bo/c, but not bo itself (see p.21jl ). 

Now, instead of u, the following 2-component spinor field $(a;) € C 2 should be introduced as the quasi-relativistic 
wavefunction of the electron: 

*(* )s (l + ^)„. (3.13) 
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In this case, integration by parts shows the time evolution to preserve the norm 



l$li 



1 + 



A 2 



A 2 



8m 2 c 2 



d 3 r v) 1 



A 2 



8m 2 c 2 



= / d 3 r u ] ( 1 + 



A 2 



d 3 r < u'u + 



At 



At 



4m 2 c 2 
0(l/c 3 ) = 



u + 0(l/c 3 ) 



2m e c / \ 2m e c 
d 3 r {u^u + Jv) + 0(l/c 3 ) = ||t/>|| 2 + 0(l/c 3 ) = 1 + 0(l/c 3 ), 



(3.14) 



while j u' f ud 3 r varies with time by a 0(l/c 2 )-amount. However, the transformation (13.13|) leaves "probability distri- 
bution" <f>^$ different from tp't/j by a fully-divergent term of the order 1/c 2 : 

= ^t^, + div )zb + 0(l/c 3 ), 
1 







ZB)i 



8m 2 c 2 



h 2 V i(u/ u) — 2ti€ijkU ! <jjP) 



(3.15) 
(3.16) 



This situation reflects the presence of negative-energy states resulting in the Zitterbewegung of the electron. 
Now expressing u and v in terms of $ using (|3.10p and (|3.13|) . write the first line of (|3 . 5|) : 



= Xu + cAv = < A + 



2m e 



-A 1 



A 



2m e c 2 



A 1 



A 2 



$ + 0(l/c 3 ). 



(3.17) 



To obtain an equation in the form ihd^/dt = h$>, one must make iterations to leave only one time derivative of $ in 
the right side of (|3.17|) . The corresponding operator is implicitly contained in A. First consider the above equation in 
the 1 /c-approximation: 



A$ = --^$ + 0(l/c 2 ). 
2m e 

After some transformations with the use of (|3.18p . we obtain in the l/c 2 -approximation: 



A + 



A 2 

2m P 



1 - 



A 2 



1 



4m 2 c 2 / 8m 2 c 2 



A, A 



A 



$ = 0(l/c 3 ). 



(3.18) 



(3.19) 



Note that commutator 



A, A 



does not contain d/dt operator, so there is only one time derivative of $ in (|3.19p . 

namely the one contained in A$. Converted into its usual form, (|3.19[) gives the quasi-relativistic equation for a 
positive-energy electron 

<9$ 



= h<5>, = * f *-divj ZB + 0(l/c 3 ); 



(3.20) 



h = 



n 



12 



2m e 
efr 



1 - 



4m 2 c 2 



2m e c 



ctH + erb + eA Q - 



r - n eh 2 tr[P[bP\] 



n 

n' 2 



4m 2 c 2 

P-b t CT, 

II 2 - 26 2 = P 2 + b 2 - 2b t a ■ P. 



(3.21) 

(3.22) 
(3.23) 



The Hamiltonian h is precisely hermitian and the corresponding equations of motion demonstrate their exact local 
gauge invariance: 

tf[A M ] = h[A^, (3.24) 
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In the 1/c-approximation, we arrive at the Pauli equation, through which the expressions for the probability current 
and density are easily found: 

<9$p 

ih ~df = /lp$P ' (3 ' 26) 
n 2 b 2 eh 

hp = (tH + eA + erb, (3.27) 



2m„ m„ 2m„c 



f p = ( C $^$ P , J- U P {P<5> P ) + (P$ P )t$p) - _^$ p(T$p | t 



(3.28) 



i.e. the current acquires an additional spin-dependent term in the 6 7^ case. The terms in (I3.27[) involving external 
fields form the interaction Hamiltonian. For A^(x) taken in the Coulomb gauge, it reads as follows: 

hpint = —A-ir + eAo- -^(t-H+-^A 2 , (3.29) 

m e c 2m e c 2m e c z 

7T = p-b t a: (3.30) 

The difference of Q3.29P from that in the conventional QED is generated by a gauge-like shift in the momentum space 
(p — ► 7r). This feature will be used in section HVl for constructing the solutions of the eigenstate problem. 

The results obtained agree with those published in [25| and [2q . in the corresponding particular cases. In the 
former paper, a nonrelativistic Hamiltonian for a free electron was obtained using the Foldy-Wouthysen method 
(l/m e -series), within the first order approximation with respect to all possible SME-corrections in the fermion sector 
of extended QED In the special case of the axial vector background 6 M , the resulting nonrelativistic Hamiltonian 

can be obtained from the formulas of paper [25j : 

hFW= 2^ e +(Tb - + 2m^^ Pl - hm) + 2r^ P {(TP) - (3 - 31) 

On the other hand, for a free electron, the Hamiltonian (|3.21[) takes the form: 

h = l^fillR ( ! _ p 2 -^-v + bi\ ah _ « 

2m e I 4m 2 ! c 2 J 2m 2 e c 2 

One can easily find that, within the linear order in 6 M and the third order in p/m e (the approximation used in 25]), 
the two expressions are identical. The absence of the term proportional to p 4 /m\ in the former expression does not 
indicate an error. Instead, it is a consequence of the difference in the expansion parameters chosen, i.e. p/m e and 
1/c, respectively. 

It should be emphasized that the method used in paper [25[ to obtain expression (|3. 3 1 [) was based on a series 
expansion of a precise relativistic Hamiltonian for a 2- component wavefunction of a free particle constructed using 
the Foldy-Wouthysen iterations 32]. Making these iterations, however, is inconvenient in the presence of external 
fields. In contrast, the method used in our paper takes these fields into account from the beginning. 

Quasi-relativistic methods similar to those used in our paper were employed in [26| to find the 1/c-corrections to 
the Dirac equation in an external electromagnetic field with additional and b^ SME-couplings. In addition, plane 
wave solutions were obtained, and SME-specific modifications of the hydrogen spectrum were estimated, within the 
nonrelativistic approximation. For the b^ coupling, the calculations performed have led to the Pauli Hamiltonian of 
the form (|3.27|) . The contributions in the fermion Lagrangian and Hamiltonian, corresponding to a M coupling, are as 
follows 0, 0]: 

A£(«) = -^o^, (3.33) 
AH {a) = 7°7 M a M =a -a-a, (3.34) 

where a M is a constant background 4- vector, which can be treated as a vacuum expectation value of some Planck-scale 
fundamental fields. As mentioned in 0, [26j], transition from the a M = to the a M ^ case is a kind of a gauge 
transformation because 

A^x) -» A^{x) = A^x) + = A^x) - d„a(x), (3.35) 
a(x) = —a^. (3.36) 
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This feature makes it possible to find a system of exact solutions of the Dirac equation modified with the a^-term 
making an inverse gauge (and unitary) transformation. Suppose the eigenstate problem is solved in the a M = case 
so that 

H^=E^\ (3.37) 
ftj 0) 4 0) =n^i 0) , £ = l,2,...,iV; (3.38) 

(\e,V£ 0) ) - / d\ ^(r)^(r) = 6 m , n , (3.39) 

where m — {mi} = {mi, mi, • ■ ■ , ttin}, n = {n{\ = {ni, U2, . . . , iif} denote the sets of quantum numbers correspond- 
ing to N hermitian operators nf^ that should commute with and with each other. The operators h\ are needed 
only to represent quantum numbers, i.e. they form a complete set of observables. It should be pointed out that the 
choice of these operators does not affect the eigenstate problem itself, but only forms the basis of the eigenstates 
and enumerates them. For example, for the nonrelativistic hydrogen atom the quantum numbers are usually taken 

such that rix = n, ri2 = I, = m define the eigenvalues of the three operators, namely h^ = H^°\ fffi = I , and 

= f 3 , in the eigenstate ipniL- 
The system of solutions for a p ^ reads 

ffW^W(r) = sC«VW(r), (3.40) 
ft< a tyW=n^W, (3.41) 

(V^,^) =Sm,n, (3.42) 

H {a) =ca(p- - (eA + a)^j + m e /3 + (eA + a ); (3.43) 

^(<0 = e la - r/hc ^\ (3.44) 
= + a . (3.45) 

The energy spectrum is shifted by a constant value ao; no spectroscopic signature is therefore left by the presence 
of the nonzero a M (i.e. transition frequencies are unaffected). However, the meaning of the quantum numbers rij 

(which run through the same set of values as in (|3.38p ) is changed, because, for tp„ functions, they correspond to 
the operators n 2 - a ' ^ which can be readily constructed from : 

For instance, in the case of a hydrogen atom, when the unitary transformation is made, if) } — > ip ) , with a being 

b^ " 

nir< 



some parameter of the transformation, the resulting V'ihL is an eigenstate of the transformed Hamiltonian H^ a \ with 

the same eigenvalue (energy), but now quantum numbers n, I, m correspond to new operators ni ^ nf \ i = 1,2,3. 
Using (|3.13[) and (|3.4[) . one can find that the transformation (|3.44j) maintains its form for the Pauli wavefunction 

$W ( r ) = e iar $ P , n (r), (3.47) 

while the nonrelativistic spectrum and hi operators arc still transformed following (|3.45p and (|3.46p . respectively. 

Another investigation, which is worth mentioning, was held in (27| . The authors have considered two non- minimal 
Lorentz- violating couplings in the fermion sector of QED: 

A£Cff.&) = ^- gv » + 5a < 75 ) 7 ^F;>, (3.48) 

where g and g a are the coupling constants while if and v% are fixed background vectors ('a' is not a component index 
but means 'axial'), and F* v = ^e^apF" 13 is the dual field tensor. We leave aside the question of the origin of such 



couplings; for more information, the reader is referred to [27J and the references therein. The authors of the paper 



also used the nonrelativistic Pauli approach to obtain the Pauli equation in the presence of the background vectors v v 
and v£, and then calculated the first-order energy corrections using perturbation theory. In addition, they considered 
an atom in a homogeneous external magnetic field also treated perturbatively. The case of the g and g a couplings, 
in general, does not permit such an easy construction of the eigenfunctions as the case of the a M coupling. However, 
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in one special case not considered in [27]], namely for the constant homogeneous external field F^ v {x) = const and 
g a = 0, the eigenstate problem can be solved in much the same way as described in (|3.44[) and (|3.45[) . Indeed, the 
transformation analogous to (|3.44|) reads as follows: 



(3.49) 



Eg) = E^+vH. (3.50) 

Again, the spectrum is shifted by a constant value, though depending on the direction of the magnetic field. However, 
the change in the wavefunctions could possibly affect, for instance, the properties of synchrotron radiation in a 
homogeneous magnetic field. 

Nonetheless, we confine ourselves to demonstrating the prospects of unitary transformations for solving wave equa- 
tions containing Lorentz-violating terms. In the following sections, similar techniques will be used to obtain the 
solutions in the case of the coupling, and to study the dynamics of a bound electron in such a background. 

IV. HYDROGEN-LIKE ATOM. QUASI-RELATIVISTIC APPROACH 

Consider first the Pauli Hamiltonian (|3.27p within the first order in b^: 

- 2 

11 eA - -^—aH + ab. (4.1) 



2m e 2m e c 
We suppose that A fl (x) is taken in the Coulomb gauge with 

divA = 0. (4.3) 
Make an inverse gauge-like shift of the momentum (II — * P) performing a unitary transformation: 

$ P <$>'p = Up<S> P , hp -► hp = Uphptfp, Up = exp {-^V ' r J> ( 44 ) 

P f eh \ 

tip = - + eA - a + fi A ) H + er ■ b, (4.5) 

Zm e \Zm e c ) 

A a = — H- (4.6) 

m e c 

It is clear that the transformation reduces the Lorentz-violating interaction to a modification of the electron magnetic 
moment, which acquires a CPT-odd correction (i A . Consequently, the terms of the first order in 6 t vanish in the 
transformed Hamiltonian as the external magnetic field H is turned off. In particular, within the approximation used, 
the eigenstate problem in an electric field would look quite conventional after the transformation. In a relativistic 
theory discussed in section [Vj an electric dipole moment correction also arises but it vanishes in the nonrelativistic 
approximation. 

Let A^ — {<p(r),0} and b — {0,0, b z }, <fi(r) being the potential of the nucleus initially considered as spherically- 
symmetric, but not mandatory the Coulomb potential. 1 The problem resembles that of an electron in a homogeneous 
magnetic field Hb but for the only difference: now there is no coupling in the kinetic term (that is, we have p for 
the momentum instead of p — -A^, where Hi, = rot A;,). The coupling to the external field b involves only the 
spin degrees of freedom but not the orbital ones. The energy eigenstates can be easily obtained in the transformed 



1 Indeed, due to radiative corrections, Coulomb attraction becomes stronger than ~ l/r at short distances contributing in the Lamb shift 
of electron eigenstates 34] . The Lamb shift however originates from the three one-loop corrections to electrodynamics including electron 
mass renormalization, its anomalous magnetic moment and the modification of the Coulomb law. For s-states, the first of them makes 
a major contribution to the Lamb shift. 
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representation and then transformed back to the initial one: 



(&p)nhn,m.(r) = Rni(r)Y^ mi (r I 'r)xm a ] (4.7) 
(*p)nta»,m,(r) = Rni{r)Y^ mi {r/r)^l+ l -^(T-r^Xrn s ; (4.8) 
E n i ms = E$+2b z m s , (4.9) 



where n = 1, 2, 3, . . ., I = 0,n — 1, m; = —I, I, m s — ±1/2 are the quantum numbers denoted according to a common 

convention; Xm a are the spin ^-component eigenvectors. R n i(r) and are the radial wavefunction and the energy 
in the b = case, respectively. In the Coulomb case, we have [32 



" ^fw^^w- <«»» 



' B 



where p = 2Zr/nrB, r-& — h 2 /m e e 2 is the Bohr radius, R = m e e A /2h 3 is the Rydberg constant, and denote the 
generalized Laguerre polynomials: 

L ( kHp) = iiP~ Ve "j^ (P v+n e- p ) , Re^>0, k = 0,1,2,... (4.12) 

The solution obtained shows that, with respect to the transformed representation, the only effect the presence of b 
generates in the leading order is a removed degeneracy over spin quantum number m s , with the energy splitting being 
< 10~ 4 Hz. Neither the spectrum nor the eigenfunctions are affected by bo, only the interaction with the external 
magnetic field is. 

The b-induced energy splitting into a doublet is a formal result of solving the eigenstate problem in a 1/c- 
approximation that does not hold true when the spin-orbit interaction is considered that removes the degeneracy 
over quantum number j. The correct splitting magnitude can be estimated by means of a perturbation theory. In 
the absence of b, the spectrum remains degenerate over I and rrij. The action of the perturbation term cr • b, however, 
preserves these quantum numbers, so the perturbation theory can be applied to the atom as to a non-degenerate 
system. This situation is typical for the anomalous Zeeman effect f3lL [33| . 

For the crb term, the energy correction was first estimated in [26| but we shall do it once again. First, let b^ = 
{0, 0,0, b z }. Following the arguments explained in the preceding paragraph, take \nljrrij) for the eigenstates in the 
b = case. Using the general expressions for them [30j, Til] , 

(r\nlj mj ) = R nlj {r)Y l jm ,{r/r), (4.13) 
Yf = 1 ( y/l + ^ + Kmj \ 

x = (_i)('- ; '+ 1 )/ 2 = ±i f or j = i± i/2 ; l' = 2j-l, (4.15) 

we obtain: 

oo 

4W r2 * • 2TTT {{l + 1/2 + Kmj) + 1/2 " xmj)) = 



oo 


that is, twice the result obtained in [261 ]. Since the corrections induced by b are minuscule, we will further treat & M as 
a purely timelike 4- vector, with the time component 6o- 

In search for ^-corrections to the eigenstates, we shall resort to the l/c 2 -approximation in the eigenstate problem. 
Consider the Coulomb case with e<j>(r) = —Ze 2 /r within the first-order approximation in 6 M = {cb t , 0}. A spectacular 
feature of this case is that the solutions can be explicitly expressed via their conventional form (for bo — 0). The 
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correspondence is generated again with a unitary transformation: 



As a result we obtain 



h = m\^oU, U = cxp {-f (l + ^) a ■ r) , (4.17) 

h - —(l- ^ \ - — + Ze2h2 (— + 2n5(r)\ (4 18) 

2m e \ 4r7i 2 c 2 / r 4to 2 c 2 \ r 3 / 



$ny m , (r) = R nlj (r) {Yj mj (r/r) - (l + (r/r)} , (4.19) 

Z 2 /ii? 



' j + 1/2 ~ An 



(4.20) 



where ;x and I' are defined in (|4.15[) . The radial functions R n ij{r) remain the same as in the bo — case (see (|4. 1 3|) ^1 . 
In the nonrelativistic limit, they take the form (I4.10|l . 

Thus, no corrections to the energy spectrum are present due to bo, within the l/c 2 -approximation. Further analysis 
will show that there are no corrections of the first order in bo (see section [Vj) . Nevertheless, the perturbative method 
used in [26[ to retrieve the energy corrections due to the term — ^?-<t • P is incorrect. The spectrum is degenerate 
over I if bo = 0, while the perturbation operator is P-odd, and hence changes the I quantum number. The expectation 
value of such an operator clearly vanishes in a state possessing a definite I, and thus definite parity P — (—1)'. The 
energy shift may not vanish, however, for some superposition of the states with opposite parities. This is common for 
the linear Stark effect [HI, [33[ that occurs due to a degeneracy of the hydrogen spectrum. The perturbation theory 
for a degenerate system must be employed instead of a simple averaging. Despite the above remarks, the methods 
employed in [2(|, have led to correct results. 

In conclusion, we will show how the correction to the magnetic moment may cause an appearance of a nonzero 
anapole moment of the atomic orbital (la |. This is a classical quantity ascribed to a parity-nonconserving system 
(such systems exist in the conventional Standard Model due to weak interaction [28| ) adding an interaction term of the 
form — Tz • rot if to the Hamiltonian of the system, with Tz being the anapole moment. Consider a hydrogen atom 
in the ground state lsi/2, mj where lower indices indicate the electron total angular momentum and its z-projection. 
Averaging the CPT-odd term — fi A ■ H in this state yields: 

Vz = <-Aa • H(r)) = - {fr A (H(0) + (r ■ V)H(O) + ...)), (4.21) 

where r = points to the center of the Coulomb field. The ground Is state possesses a definite parity (—1)' = +1 
and, in addition, a spherical symmetry, hence 

(Aa> = 0, (4.22) 
(xiX k ) = ^S ik {r 2 ), i,k = 1,2,3. (4.23) 

With the help of expressions (|4.7ll and (|4.10p . one can easily find that 

(aiX k x n ) =2r^d i3 8 kn ■ nij, i, k, n = 1, 2, 3. (4.24) 

which results in the following: 

V z » -^fm J e 3lk d l H k = -Tz-rotH, (4.25) 

m e c 

Tz = 2erl (-^) m 3 e 3 , (4.26) 



where e 3 is the basis unit vector along the z-axis. 



V. SERIES EXPANSION OF THE DIRAC EQUATION WITH RESPECT TO b 

In this section we discuss the case b» = {b°,0} and = {A ( c) (x) + <j>[r), A {c \x)} using the Heaviside units, 
with h = c = 1, a = j-. Consider the Hamiltonian (12. 5|) and transform the corresponding wave equation using the 



10 



gauge-invariant unitary transformation: 



ip(x) -> tp(x) = e-^^ipix), (5.1) 

^-4-^-4= e "^(^-4) e ^ (5 - 2) 

A A = S.f- — (S-i + l) 7 °75, (5.3) 
m e 

£ = [rP] = -[Pr]. (5.4) 
Restricting ourselves to the second-order approximation in bo, we obtain: 

H w a(p — ei< c ») + /3m e + e(0 + 4 c) ) - —/V - ^-E (c) - Aa-^ (c) + S^[A^], (5.5) 

with / = + 1. ff^f^] stands for the second-order terms in bo describing the interaction with the external field 
. Additional electric and magnetic dipole moment operators read as follows: 

= — 7°PH (5.6) 

m e 

= -i75AA = -— -W- (5-7) 

As we can see, no non-linear terms in the external field are present up to the first order in bo, inclusively. Moreover, 
the moment (1a couples with the external field only but not with the spherically-symmetric 'background' field <p(r), 
because for such a field (La ■ (— V0) = 0. The same situation holds in every higher order of the expansion due to the 
fact that [A^,<^(r)] = 0. For the same reason, the expressions for operators cLa and fi A are n °t affected by 4>(r), in 
particular, they maintain their form for a free electron. Our approach however is applicable only to systems with the 
effective size much less than 1/bo > 10 _3 cm. 

We did not obtain a CPT-odd correction coupling to the electric field in the 1/c-approximation because, in contrast 
to fi A , is a block-off-diagonal operator, 

d A = ( r ° 1 [ °? ] ) , (5.8) 

m e c \~[<tx] J v > 

which mixes the 'upper' and the 'lower' 2-component spinors of the wavefunction. The 'lower' spinor vanishes in the 
non-relativistic limit (see eq. (|3.10|) ). and so does the operator cLa- Instead, ft, A is a block-diagonal matrix that mixes 
the 'upper' spinors with themselves and consequently it does not vanish in the nonrelativistic limit. 
Let — 0, then the spherical symmetry allows us to search for the eigenfunctions in the form 



* , (rt)-( <V r) ^ (r/r)> 



(5.9) 



where n r = n — j — 1/2 is the radial quantum number and n is the principal quantum number, and I, as usual in the 
relativistic theory, determines the parity of the state P = (—1)', but not its orbital momentum. Operators (1a and 
fx a have vanishing expectation values in such a state, ip is the eigenfunction of the operator /7 , with the eigenvalue 
/ ee x(j + 1/2). In the case under consideration, the transformed Hamiltonian (|5.5p is the sum of its conventional 
value (for bo — 0) and a term proportional to /7 ; ip is an eigenfunction for both of them if the radial functions R( u ' v > 
are taken the same as those for the bo = case. The energy value which responds to ip is 



E = E = E% - x(j + 1/2)-^ = E^. ± (j + 1/2)-^ for l = j± 1/2. (5.10) 

r J m e ' J m e 

An additional foo-induced second-order energy splitting therefore arises: 



b 2 

AE(j) ee E nrJ+1/2J - E nrJ _ 1/2d = (2j + l)-2-. (5.11) 

7Tl e 
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This term originates from parity violation due to the 6^-induced violation of CPT and removes the degeneracy over 
I in the Coulomb field case. \AE(j)\ < 10 5 Hz for j = 1/2, that is, four orders of magnitude smaller than the Lamb 
shift [34]. Nonetheless, in contrast to the latter one existing mainly for s-states, the splitting (|5.11|) increases with 
growing j. 

The eigenfunctions in the initial representation are obtained after performing the inverse transformation: 

iiniim(r)=e- b ° f / 2m *e- b ° r / 2 \ m \ m ' I m . (5.12) 

The presence of the admixture of spherical spinors with the different value of the orbital quantum number (I') breaks 
the parity of the states. 'Probability distribution' i/j'tp is not affected however, compared with the conventional 
(bp = 0) case, within the chosen approximation. Since the conventional solution in the Coulomb case is well-known 
[30(, we can explicitly find the second-order approximation for the eigenfunctions in the bo ^ Coulomb case: 

<\ =- e (l+f 7 4M I"' 2 , (5.13) 

M )n - \ 1/2 

n r j) n r- \ Xr 




R(u) X = ( (me±Qv! x 1 1 

R {v) J 1 ' 4m e ^(/ + ^)r(2 7 + TV + l 



*(2Ar)" 



x ((f + ^) 4 2 ; } (2Ar) ± (1 - 5„ r , )(2 7 + n r )4 2 ;2 1 (2Ar)) , (5.14) 



A = V m e-C 7-\/(.7 + l/2) 2 -(Za)2, (5.15) 

with being the generalized Laguerre polynomials defined in (|4.12p .The expression for the energy demonstrates an 
alternative mechanism of removing the degeneracy over I, different from that connected with the one-loop corrections 
in quantum electrodynamics. 



VI. SPECIFIC RADIATIVE PROPERTIES OF A HYDROGEN ATOM INDUCED BY b 



Finally, we demonstrate an example with CPT and Lorentz violation leading to radiative effects specific for the 
bo case and linear in bo- Following the system of units convention used in section [TTTl we assume b^ = {c6 t ,0}. 
Since the primary goal of this section is to obtain the leading-order 6o-induced terms in the radiation distribution, 
we restrict ourselves to the Pauli approximation and consider the radiation of a hydrogen (Z = 1) atom. Upon the 
transformation (|4.5[) . the only term remaining with bo is —fij^ ■ H. This term violates the spatial parity of the atom. 
With the use of the standard formulas [3l| . we find the angular distribution of spontaneous radiation probability: 



dW fi (k,X) 



\k\ 
ill 



uj/c = (Ei -E f )/Hc>0, 



er — — (k ■ r)r — 
(I + a) + p, A , 



k 



x /x 



2m e c 



A = l,2; 



(6.1) 
(6.2) 
(6.3) 

(6.4) 



where k, A are the photon momentum and polarization, and e^(k) is the polarization vector. (/| and \i) denote the 
final and the initial electron states. The correction fi A to the magnetic moment operator fi is defined in (|4.6p . 

Radiation processes allowed due to the parity-violating interaction with fi A (further referred to as Aj-radiation 
with j denoting the photon angular moment) are restricted by the same selection rules as those for the _El-radiation, 
the corresponding matrix element having the form typical for Ml-radiation. Thus Al- and .El-photons have the same 
multipolity but the opposite parity. Linear in bo corrections to the angular distribution occur due to the interference 
between the El- and the Al-radiation. Consequently, they vanish over the whole sphere because of the cancellation 
of spherical spinors with different parities, resulting in no linear in bo terms present in the total transition rate. 
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FIG. 1: Angular distribution of spontaneous radiation for 2pi/2,i/2 — ► lsi/2,-1/2 transition 



The interference term does not vanish, for example, for the transition 2p 1 / 2 ,i/2 — > l s i/2,-i/2- The calculations show 
that after averaging over the photon polarizations A = 1,2, the resulting angular distribution reads (0 is the angle 
between k and the z-axis): 



l + cos 2 9 + \ cosG \. (6.5) 



dW _ 512a 3 R 
dflk 6561-7T 

As we can see, the presence of bo induces the violation of the conventional 'fc-parity' of the distribution (the radiation 
rates in the opposite directions differ in the bo 7^ case). The relative magnitude of this violation is of the order 
\bo\/m e c 2 < 2 • 10~ 8 . Distribution (|6.5p is depicted in figfTJ with the dotted curve related to the bo = case. To make 
the picture more vivid, we chose bo/rn e c 2 = 0.05. 

For unpolarized atoms, i.e. after averaging over rtij, rrij' quantum numbers, the spherical symmetry is restored 
in the distribution, with no linear in 60 k-odd contributions present. This is the consequence of S'0(3)-invariance 
unbroken even in the bo 7^ case (while 0(3) symmetry is broken since 60 is a pseudoscalar). 

We left aside the problem of polarization of atoms. If one uses Zeeman effect in a homogeneous magnetic field to 
obtain the polarization, then this magnetic field would also lead to parity violation due to the interaction with fi A . 
Another way the external magnetic field can break the atomic P-parity is that in the reference frame of a moving 
atom, an additional electric field will be induced that breaks the parity. The distribution of radiation of moving atoms 
can also be shifted due to aberration. These problems need further consideration. In the present paper however, we 
just demonstrated yet another scenario of P-parity violation in atomic transitions. 



VII. CONCLUSION 



In this paper, we have considered several solutions of the Dirac equation in the framework of the Standard Model 
Extension with particular types of Lorentz violation. The l/c 2 -approximation for the extended Dirac equation was 
derived in the background of the axial vector SME-coupling & M . The expansion of the relativistic Dirac equation with 
respect to bo has been employed to solve the eigenstate problem for an electron in a spherically-symmetric potential 
well. The unitary transformation was found that was used to express the solutions with bo 7^ in terms of solutions 
for bo = 0, with the second order accuracy with respect to bo- Explicit solutions have been obtained in the case of the 
Coulomb potential, demonstrating a specific 60-quadratic energy splitting. The degeneracy over the orbital quantum 
number is removed, and it was shown that the corresponding energy splitting does not vanish for large j. 

In addition, unitary transformations were used to obtain the exact eigenstates in the case of the coupling 
—g-tpj^v^F*^ with constant homogeneous electromagnetic field F^ v p^ |. 

The unitary transformation made it possible to obtain the ^-corrections to the operators of El and Ml moments 
of the electron. These moments effectively lead to an existence of the anapole moment of the orbital [HI . 

Finally, the distribution of spontaneous radiation of a polarized hydrogen atom was shown to lose its central 
symmetry in the bo 7^ case, due to the violation of spatial parity. The results obtained can be treated only as an 
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illustration of the application of the model adopted. There are other physical effects that should also be considered 
together with the one discussed in this paper. 
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